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Abstract We study the refractive-index sensing properties of plasmonic nan- 
otubes with a dielectric core and ultra-thin metal shell. The few-nm thin metal 
shell is described by both the usual Drude model and the nonlocal hydrodynamic 
model to investigate the effects of nonlocality. We derive an analytical expression 
for the extinction cross section and show how sensing of the refractive index of 

ryj ■ the surrounding medium and the figure-of-merit are affected by the shape and 

size of the nanotubes. Comparison with other localized surface plasmon resonance 
sensors reveals that the nanotube exhibits superior sensitivity and comparable 

2~H' figure-of-merit. 
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1 Introduction 



It is well known that metallic nanoparticles can sustain localized surface plasmon 
(LSP) oscillations, whose resonance frequencies in the quasi-static limit depend 
solely on the geometry of the nanoparticle, the permittivity of the metal and the 
surrounding permittivity. The dependency of the LSP resonance (LSPR) on the 
l^-s ' surrounding medium makes metallic particles extremely good sensors, progressing 
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towards the detection of single molecules 1 . However, the weak effect of retarda- 
tion on the LSP resonance in nanosized metal particles leaves only one parameter 
to truly engineer: the geometry. By modifying the structure of the metal nanopar- 
ticle to have a dielectric core with a metal shell, an increased tunability is achieved 
due to the plasmon hybridization of the inner and outer surfaces of the metal [2] . 
Especially the spherical core-shell structure has received a considerable amount of 
attention in recent years [HHHEM] due to its excellent and tunable sensing proper- 
ties, which show great promise in biological studies such as cancer therapy [7|. The 
plasmon hybridization allows one to position the LSP resonance of the nanoshell as 
desired by simply varying the core size n and/or outer radius ri appropriately [8] . 

The hybridization of the inner and outer surface plasmons increases when the 
metal shell becomes thinner [8], which gives rise to significantly altered LSP res- 
onances compared to usual homogeneous metal nanoparticles. Studies of the hy- 
bridization between two spherical [9] or cylindrical [10] metal nanoparticles in 
few-nm proximity reveal that effects of nonlocal response increase with increas- 
ing hybridization. Furthermore, nanosized metal particles |HU12|[T3lll4] and metal 
films [15 are also strongly affected by nonlocal effects. The core-shell particle thus 
calls for a nonlocal description, since it features an ultra-thin metallic shell with 
resulting strong plasmon hybridization. 

The use of arrays of nanotubes with high aspect ratio for biosensing [16] and 
hydrogen sensing [17] has yielded impressive results, yet only few theoretical stud- 
ies have been performed on the nanotube [181119] . Schroter et al. investigate the 
plasmonic modes and dispersion relations of the nanotube [18 , while Zhu et al. 
perform calculations using the discrete dipole approximation to discuss the changes 
of the resonance wavelength of the nanotube due to variations of the aspect ratio 
[19] , Thus, to our knowledge no systematic study has yet been performed that 
addresses which parameters determine the LSPR refractive-index sensitivity of a 
nanotube-based sensor. In this paper, we fill this gap with a systematic study of 
the sensing and scattering properties of a single infinitely long cylindrical core- 
shell nanowire (see inset of Fig.rj}, which is a good description of dilute arrays of 
non-interacting nanotubes with high aspect ratio. On the basis of this study, we 
propose how to optimize a nanotube-based sensor to achieve the utmost sensitivity 
for the refractive-index sensing of both gases and liquids. 

The outline of this paper is as follows. In Sec. [2] we discuss the physical prin- 
ciples of local and nonlocal response, and introduce the sensitivity and figure-of- 
merit (FOM) as quantitative measures of the performance of a LSPR-based sensor. 
Section[3]is dedicated to the study of a nanotube with a silica core and gold shell. 
We determine the dependency of the sensitivity and FOM on the shape and size 
of the nanotube, using both local and nonlocal theory to model the response of 
the gold shell. Our conclusions and outlook on nanotube-based sensors is given in 
Sec. [31 and details on the analytical calculations in the Appendix. 



2 Theory 

The ability of LSPR-based sensors to detect changes in the refractive index of their 
surrounding medium is usually quantified by the sensitivity and FOM [1 . The 
sensitivity dX/dn^, is determined as the shift in wavelength of the considered LSP 
resonance in the extinction spectrum of the sensor, when varying the background 
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refractive index rib = y/^b, while the FOM is given as 



F0M= M9n^ 

AX y ' 

where AX is the resonance linewidth, calculated as the FWHM of the considered 
LSP resonance in the extinction spectrum. Thus, to determine the performance of 
the nanotube as a LSPR sensor, we must calculate its extinction cross section, as 
this quantifies the extinction spectrum and therefore allows us to determine the 
sensitivity and FOM. 

Predictions for the extinction cross section depend on how the optical response 
of electrons in the metal is modeled. The common approach to describe the re- 
sponse of metals is by making the local approximation which assumes that the 
response field at a certain position is proportional to the driving field at that posi- 
tion, with the proportionality function being a position- and frequency-dependent 
dielectric function. This approach has the rather unphysical consequence that all 
surface charges reside on an infinitely thin layer on the boundaries of the metal, 
thereby neglecting the actual extent (or wave nature) of the electrons. While the 
local approximation is justified as long as the metal boundaries are far apart such 
that the interaction between electrons due to their extent can be neglected (i.e. 
large metallic structures) , it can not be safely assumed for nanosized metal parti- 
cles where the wavelength of the electron becomes comparable in size to the metal 
particle. Describing the metal using the semiclassical hydrodynamic Drude model 
|14j . we relax the local approximation by allowing the existence of local inhomo- 
geneity in the density of the electron gas, which gives rise to pressure waves. The 
electron-gas pressure waves provide a means to transport energy in the metal in 
addition to the electromagnetic waves, which gives rise to nonlocal response: the 
response of the metal at a certain spatial point can depend on the driving field at 
other nearby points (on the length scale of the Fermi wavelength) in the metal. 

In the Appendix, we provide an analytical expression for the extinction cross 
section in the cases of both nonlocal and local response, for a normally incident 
TM-polarized wave, see the inset of Fig. [TJ We have checked the analytical expres- 
sion with our numerical implementation of the hydrodynamic Drude model |10j . 
which showed perfect agreement (not shown in this paper). 



3 Results and discussion 

We consider the specific core-shell structure, where the core is silica (SiC>2) with 
dielectric constant e c = 1.5 2 and the shell is gold (Au) modeled with the data by 
Rakic et al. [20] . To clearly show the difference in extinction cross section in local 
and nonlocal response, we start by examining the case where interband effects in 
Au are neglected. Figure [1] depicts the extinction cross section for a (ri,ra) = 
(40nm,45nm) silica- Au cylinder in vacuum comparing the local and nonlocal 
model. The local approximation shows three distinct peaks, two at low frequencies 
(dipole and quadrupole peaks) and one at a high frequency (near 7eV). These 
are due to the interaction between the localized plasmons at the inner and outer 
surface of the nanoshell, or equivalently, the interaction between a cavity mode 
and a cylinder mode [8]. The nonlocal description allows the same classification 
of peaks as the local approximation [141110] , although the high-frequency peak is 
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Fig. 1 Extinction cross sections as a function of incident photon energy for TM-polarized light 
normally incident on a (n , ra) = (40 nm, 45 nm) silica- Au cylinder in vacuum. The three curves 
correspond to the nonlocal and local models without interband transitions (solid green and 
dashed blue curves, respectively), and the nonlocal model with interband transitions (solid red 
curve). Free-electron parameters for Au as in Ref. |20| : 7kj p = 7.872eV, h-f = 0.0530eV, and 
tip = 1.39 X 10 6 m/s. Interband parameters for Au are also as in Ref. 20 and valid up to 5 eV. 
The panel on the right shows the normalized intensity distributions |E| 2 /|Erj| 2 in the nonlocal 
model without interband transitions at the dipole and quadrupole resonance frequencies. Here, 
Eo is the incident electric field. Inset: Schematic diagram of core-shell structure with relevant 
parameters. 



blueshifted compared to the local model. Since sensing depends on peak shifts, it 
is important to take possible nonlocal blueshifts into account. However, the low- 
frequency resonances show no noticeable blueshift, because the strength of the 
nonlocal blueshift does not only increase with decreasing thickness of the metal 
layer 14,9, 15, 13 1 but it also depends on the frequency, with a decreasing blueshift 
for lower frequencies. Thus we find that there is an intricate interplay between 
plasmon hybridization and nonlocal response: Since a thinner metal shell produces 
stronger plasmon hybridization, the dipole and quadrupole peaks are pushed to 
such low frequencies that the nonlocal blueshift effect due to nanosized metallic 
features is counteracted by the low frequency of the resonances. 

The panel on the right of Fig. Q] shows the nonlocal normalized intensity distri- 
bution in the metal at the dipole and quadrupole resonance frequencies, illustrating 
the expected dipole and quadrupole nature of the resonances. Above the plasma 
energy ftu p we see the characteristic additional resonances in the nonlocal model 
due to the excitation of longitudinal modes, as previously reported for different 
metal nanoparticles |131I141[21] . 

The difference between the red and green curves in Fig.LTJshows the importance 
of taking into account interband transitions in the response of the metal shell. The 
implications on the dipole and quadrupole resonances are that they are redshifted 
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and damped due to interband transitions, with greatest impact on the quadrupole 
peak. In the remaining part of this paper, we will always use measured values for 
the dielectric function [20J, i.e. we take interband transitions into account. We will 
concentrate on the dipole resonance, since this peak is the strongest, is close to 
visible and infrared frequencies and can be affected by the shape and size of the 
cylinder and the background permittivity. 

There are two geometrical properties that can be modified in the nanotube 
structure: the first is the shape defined by the r\/r2 ratio and the second is the 
overall size, that is, varying the outer radius r-2 but keeping r±/r2 constant. In Fig. 
[2] we show the effect of shape variations of the nanotube on its sensing abilities, 
which is quantified through the change in the dipole resonance wavelength when 
the background refractive index is increased. We see that regardless of the shape, 
the dependency is always approximately linear. However, as shown in Fig. [2ji) 
there is no significant dependency on the background refractive index for low 
ri/r2 ratios, indicating the lack of ability to sense. Only when the shell becomes 
thin (n/ra —> 1) does the resonance wavelength shift with the refractive index. 
The thinner the shell, the greater is the average slope of the curves. Relaxing 
the nonlocal description to a local one does not change this trend, because the 
dipole resonances occur at too low energies for the nonlocal blueshift to kick in. 
Furthermore, the resonance wavelength shifts to higher wavelengths when the shell 
becomes thinner, because the coupling between the cavity and cylinder modes 
increases. Thus, even though Fig. [2jiv) represents a nanotube with a 2 nm thin 
metal shell, where nonlocal blueshifts are expected to be very prominent, the 
local approximation predicts sensitivities that are almost identical to the nonlocal 
description. So, as in Fig. [1] here in Fig. [2] we see that for ultra-thin nanotubes 
the usual observation of larger nonlocal blueshifts for smaller structures does not 
occur. The nonlocal blueshift cancels out with the decrease of the resonance energy 
due to increased hybridization. 

For a more quantitative description of the sensitivity of the nanotube, we 
present sensitivity and FOM calculations of the nanotube structures shown in 
Fig. [2] at the refractive index of air and water in Table [1] As in Fig. [2] it is 
again clear from Table [1] that increased sensitivity can be achieved for thinner 
metal shells. Comparing the sensitivity of the nanotube with other LSPR sensors 
based on different nanoparticle geometries [22] > where the sensitivity is in the range 
90—801 nm per refractive index unit (RIU), shows that the nanotube is comparable 
in sensitivity for ratios ri/r2 > 0.7, while it is superior for very high ri/r2 ratios. 
Comparison of the FOM with other nanoparticle LSP sensors also shows equally 



Table 1 Sensitivity and figure-of-merit calculations Eq. Q in the nonlocal description at the 
refractive index of air nj, = 1 (for gas sensing) and water n^ = 1.333 (for liquid sensing) for 
the four different shapes of Fig. [2] 
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Fig. 2 The dipole resonance wavelength calculated with both local and nonlocal response 
taking into account interband transitions as a function of the background refractive index for 
four different r\/r2 ratios: (i) 0.7 (ii) 0.9 (iii) 0.95 (iv) 0.98. The outer radius of the nanotubc 
is kept constant at T2 = 100 nm. The panel on the right shows the normalized intensity 
distribution |E| 2 /|Erj| 2 in the nonlocal model at the vacuum dipole resonance wavelength for 
the corresponding four different shapes. 



good performance by the nanotube, although the FOM is mainly dependent on 
the properties of Au and not easily improved by changing the geometry 23 1 . The 
sensitivity values in Table [1] also reveal that the nanotube has a high sensitivity 
at both the refractive index of air and water, which shows the versatility of a 
nanotube-based sensor and its applicability as both a gas and liquid sensor. 

Besides shape variations, we also varied the size ri of the nanotube, while 
keeping r\/r2 constant. Figure [3] depicts the dipole resonance wavelength as a 
function of the background refractive index for three different sizes with r\/r2 = 
0.9. The sensing ability of the nanotube is not as dependent on size as it is on shape, 
which can be seen by the three almost parallel lines in Fig. [3] Even though the 
sensitivity does not change much with increasing size, there is still an optimum 
size which occurs at ri = 50 nm and ri = 70 nm for liquid and gas sensing, 
respectively, see the inset of Fig. [3j The fact that it is neither the smallest nor 
the biggest nanotube size that gives the highest sensitivity can be explained by 
a trade-off between the total structure size and the shell thickness. If the size of 
the structure is too small, then we have a weak LSP excitation and thereby poor 
sensing ability, but if the structure size is too big (with the shape kept constant) 
the absolute shell thickness increases, which also decreases the sensing ability, as 
we saw in Fig. [2] Therefore, for a larger r\/r2 value the optimum size will also be 
larger. 

In Fig. [3] we also show the calculations using the local approximation. As seen, 
effects are surprisingly well accounted for even with a local description, despite 
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Fig. 3 The dipole resonance wavelength calculated with both local and nonlocal response 
taking into account interband transitions as a function of the background refractive index 
for three different r% values: (i) 30 nm (ii) 50 nm (iii) 80 nm. The shape of the nanotube is 
kept constant by setting r\/r2 = 0.9. The panel on the right shows the normalized intensity 
distribution |E| 2 /|Erj| 2 in the nonlocal model at the vacuum resonance wavelength for the 
corresponding three different sizes. Inset: The LSPR sensitivity at the refractive index of air 
(rib = 1) and water (nj, = 1.333) calculated with the nonlocal model as a function of outer 
radius while keeping r\/T2 = 0.9. 



the fact that we actually consider very thin metallic shells, for instance a 3nm 
shell in Fig. [3](iii) , with concomitant strong plasmon hybridization. The strong 
hybridization in ultra-thin metal shells shifts the dipole resonance to very low 
energies, where the nonlocal blueshift is weak. The sensitivity and consequently 
the FOM are therefore weakly influenced by nonlocal response. Although it is 
hardly visible in Fig. [3l the local resonances do in fact occur at slightly longer 
wavelengths than in the nonlocal description, revealing a small nonlocal blueshift. 



4 Conclusions and outlook 



We have examined the infinite single dielectric-metal nanotube structure as an ap- 
proximation for a dilute array of nanotubes with high aspect ratio. We calculate 
the extinction properties of a silica-gold nanotube analytically for both local and 
nonlocal response by extending the Mie theory for nanowires to nanotube geome- 
tries. Our investigation reveals that in contrast to the spherical nanoshell [6], the 
sensing ability of the nanotube is highly dependent on the shape of the structure, 
where few-nm thin shells produce extreme sensitivities. The sensitivity is shown 
to be less dependent on the overall structure size. The sensitivity at the refractive 
index of air and water of ultra-thin nanotubes are superior to other nanoparticle 
geometries, making nanotubes very promising for both gas and liquid sensing. 
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Our results also show unexpectedly that nonlocal response has negligible in- 
fluence on the extinction and sensing properties of the nanotube, even though the 
metal shell is ultra-thin (a few nm), because the hybridization in the nanotube is 
so strong that the dipole resonance is pushed to very low energies. The strength 
of the nonlocal blueshift is an interplay between the metal thickness and the res- 
onance energy, where a thinner shell produces a stronger blueshift while a lower 
energy produces a weaker blueshift. This interplay is surprisingly well-balanced 
in the nanotube structure, because a thinner shell gives rise to lower resonance 
energies. 

With the high sensitivity and good FOM of the nanotube geometry, we propose 
a sensor based on ultra-thin nanotubes. The robustness of the sensitivity of the 
nanotube to size variations provides desirable advantages, since fluctuations in 
size due to imperfect fabrication will have less impact. In the special case of gas 
sensing, the sensitivity may be further improved by a factor of two by designing the 
nanotube to have a hollow core. With a hollow core, the inner surface of the metal 
shell is also exposed to the surrounding medium, which significantly improves the 
sensitivity. However, mechanical stability is sacrificed with a hollow core if for 
instance the nanotubes are to stand vertically on a substrate. 



5 Appendix 

The nonlocal optical properties of the nanotube are determined by solving Maxwell's 
wave equation coupled to the hydrodynamic equation for the current [14 . We solve 
the coupled set of equations by extending the Mie theory for wires of Ref. [23] to 
core-shell structures. By expanding the electromagnetic fields in the dielectric core, 
metal shell and surrounding medium in cylindrical Bessel functions, we can most 
easily take into account Maxwell's boundary conditions along with the additional 
boundary condition of a vanishing normal component of the current in the non- 
local case [T2J. Although quantum tunneling is not taken into account with this 
treatment, we do not expect any such effects to be important in this structure [251 

To determine the extinction property of the infinite cylindrical nanotube we 
calculate the extinction cross section [27] 

2 °° 
CTeXt = "fc^ ^ Re{a„}, (2) 

n= — oo 

where ko = y/e^ui/c is the background wave vector, £b is the background permit- 
tivity and a n is a cylindrical Bessel-function expansion coefficient for the scattered 
electromagnetic field. We consider a normally incident electric-field polarization 
perpendicular to the cylinder axis (TM), as sketched in the inset of Fig. [I] The 
nonlocal-response scattering coefficient is calculated analytically as 

e^J n (kor2) [C n + J n P n - H' n Q n ] - \feJ' n (kor2) [JnPn, - H n Q n ] 



~^H n {k r 2 ) [C n + JkPn - H n Q n ] - y^H n {k r 2 ) [J n P n - H n Q n ] ' 



(3) 

Here, J n and H n are the Bessel and Hankel functions of the first kind, fc t = 
y/eu/c and e(w) = e ther(^) — w p /(w[w + vy]) is the Drude local-response function 
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that includes interband effects through e ther(k>). The argument of the Bessel and 
Hankel functions are k\r 2 unless written explicitly otherwise. 
The coefficients P n , Q n and C n are given by 

P n = p n a n + J„(fc c ri) [H n (k t ri)5 n + H n r n ] , (4) 

Qn = q n a n + J n (fc c ri) [J n (hri)5 n + J n r n ] , (5) 

ZTi 

Cn = [Hn(kir 2 )c n - Jn{hr2)dn] , (6) 

K0T2 

where k c = y/e^u/c and e c is the dielectric constant of the core. Furthermore, 
fc[ = (uj 2 + iu"f — ojp/e ther)//3 2 and (3 2 = 3up/5 with vf being the Fermi velocity 
of the metal shell. The coefficients p n , q n , 0t n , 5 n and t„ of Eqs. ([3][5]) are given as 

Pn = y/eJn(k c ri)H n (ktri) - y/eZJn(k c ri)H' n (k t ri), (7) 

q n = y/eJnikcr^Jnihr!) - y/€ZJn{k c r 1 ) J n (ktri). (8) 

a-n = ( X ° ther j [Jn{ k i r 2)H' n (kiri) - H' n (kir 2 )J n (kiri)] , (9) 

< - '''"'" v ~' '"" 1 "" e °-^l[4(fe 1 r 2 ) J ff„(fc 1 r 1 )-<(fc 1 r 2 )J n (fc 1 r 1 )] J (10) 



k t k 2 r 2 

k\Tl -^/ecCotherCe — Eother) 



[H' n (hri)Mhn) - f n (kiri)H n (kin)] , (11) 



ktk 2 nr 2 
while the coefficients c n and d n of Eq. ([6]) are given as 

Cn = fn [Jn(for2)rj„ + J n (kin)K n ] + j' n {k]ri)g n [JnPn - H n q n ] , (12) 

d n = fn [H' n {kir 2 )rj n + H n (k i r 1 )K n ] + H' n (kiri)g n [JnPn - H n q n ] , (13) 

where 

mfcl£other(e - Eother) ,-,.-, 

9" = ri. ' ( 14 

koktr 2 

U = in ^\- eothei) J n (hn), (15) 

r] n = k\ [J n (hri)Hn(hri) - H n (k t r 1 )J n (k t r 1 )] , (16) 

_ n (e - Eother) 



fc t r 2 ri 



[Jn(*tri)fl-„-H„(fctri)J„]. (17) 



The local-response result can be retrieved in the limit of a vanishing Fermi velocity 
for which P n = p n , Q n = q n and C n = 0. 
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